Still water: dead zones and collimated ejecta from the impact of granular jets 
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When a dense granular jet hits a target, it forms a large dead zone and ejects a highly collimated 
conical sheet with a well-defined opening angle. Using experiments, simulations, and continuum 
modeling, we find that this opening angle is insensitive to the precise target shape and the dissipation 
mechanisms in the flow. We show that this surprising insensitivity arises because dense granular jet 
impact, though highly dissipative, is nonetheless controlled by the limit of perfect fluid flow. 
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Students are familiar with liquids as an intermediate 
state of matter: like gases, they flow easily but, like 
solids, they are condensed due to inter-particle attrac- 
tions. Later they may be taught that liquids can be 
modeled without attractions if the particle density is kept 
high by confinement [1-3 . However, even without attrac- 
tions or confinement, non-cohesive particles can behave 
like a liquid: when a high-density jet of grains hits a 
target it ejects particles in a thin sheet similar to that 
created by the impact of a liquid jet [4l l5l. Here, we 
investigate why these two different types of materials be- 
have in such a similar way. 

In our experiments we measured the velocity near the 
impact center of the dense granular jet. Our measure- 
ments reveal a dead zone, a region of nearly immobile 
particles, instead of the smoothly varying straining flow 
characteristic of Newtonian liquid impact. We then used 
discrete-particle simulations to examine how this quali- 
tative change in the velocity field alters the ejecta and 
found that it is only weakly modified. To understand 
the origin of this insensitivity we first used experiment 
and simulation to show that the granular motion in this 
regime can be modeled by an incompressible frictional 
flow, then analyzed the behavior of the ejecta in the con- 
tinuum model in the limit where the dissipation vanishes. 
We find that the form of the ejecta is dictated by inertia 
and therefore highly robust. The angle of the emerging 
ejecta is insensitive to the type of dissipation mechanisms 
present, whether the central region is flowing or static, 
and the target shape. 

This study shares with experiments at the Rela- 
tivist ic Heavy Ion Collider (RHIC) an interest relating 
the ejecta pattern to bulk properties inside the impact 
zone 0. Those experiments show that the collisions of 
gold ions produce surprisingly collimated scattering pat- 
terns. Some researchers have interpreted this coherent 
ejecta as evidence that the quark gluon plasma forms a 
dense, nearly perfect liquid [THO]. Granular jets that also 
scatter like liquids constitute a macroscopic analog. In 
granular physics, low-speed, dense flows [TOl[TT] and high- 
speed, dilute flows p!2HT5 have received much attention, 
but less is known about the high-speed, dense regime we 
examine here. This regime is relevant in a wide variety 
of contexts. Two examples are protoplanetary formation 
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FIG. 1. Experiments reveal that a large dead-zone forms dur- 
ing granular jet impact, (a) Schematic: a burst of air ejects a 
dense column of spherical glass beads from a long tube. The 
granular jet subsequently hits a target at a speed Uq. (b) Im- 
pact against the target produces an axisymmetric ejecta sheet 
with opening angle ^o- Long-time exposure of granular jet 
impact as viewed from (c), the side, and (d), below, reveals 
a large dead-zone. 



by the collision of dust aggregates [TF- 20 , and abrasive 
blasting using high speed sand jets [21^,22^. 

In our experiment, depicted in Fig.[T](a), we follow the 
protocol of Cheng et al [4|. High-pressure gas pushes 
a dense plug of non-cohesive glass beads of radius Rq = 
54±9 /im out of a tube of radius i?jet =0.8 cm at a speed 
Uq. Depending on the gas pressure, Uq is between 1 m/s 
and 16 m/s. This jet of particles hits a target of radius 
i^Tar- The particles are then ejected from the target at 
an angle ^o in a thin axisymmetric cone (Fig. flVb)). 

We present measurements of the internal flow for im- 
pact onto a roughened steel target with radius i^xar = 
i^jet and onto a smooth, transparent glass target with 
radius i^xar = 6.4i?jet. We view the impact zone near the 
steel target by slicing the originally axisymmetric exper- 
iment in half along its length and observing the central 
region from the side through a glass window (Fig.flVc)), 
and view the impact zone near the transparent target by 
observing from the rear (Fig. flVd)). The interior of the 
jet, as viewed from both the side and rear, reveal a region 
extending over a signiflcant fraction of the target where 
the grain motion is negligible compared to the surround- 
ing flow. We call this the dead zone. 
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FIG. 2. PIV measurements of the azimuthally averaged 
radial speed of grains at the transparent target {ur{r, z = 0)) 
measured from below. We see a dead zone at the target center. 
Inset shows the effective granular temperature Teff, or the 
velocity fluctuations, normalized by Uq. The dead zone is 
not only static but also cold. 



Fig. |2] shows the azimuthally- averaged radial veloc- 
ity along the smooth transparent target. Using particle 
tracking, we also measured the velocity fluctuations along 
the target base to obtain the normalized azimuthally- 
averaged granular temperature {T^^) /Uq (Fig. [2] inset). 
The normalized temperature measures how much energy 
is contained within the velocity fluctuations relative to 
the kinetic energy of the incident jet. Its small numerical 
value shows that the kinetic energy originally possessed 
by the now immobile particles in the dead zone is largely 
dissipated. Fig. [3](a) shows speed contours in front of 
the roughened steel target, where r and z are the radial 
and axial coordinates. Defining a particle as being in 
the dead zone when \\i\^ /Uq < 10~^, we see that this 
impact produces a dead zone that is broader and taller 
than the one obtained with a smooth target. Specifically 
the dead-zone radius i?DZ is 0.76i?Tar while the height 

Hy^Z is 0.4i?Tar. 

In the dilute jet limit, the ejection angle ^o of a sin- 
gle particle changes dramatically if the particle collides 
with a fiat target instead of a conical dead-zone. While 
Cheng et al. [4 first noted that dense granular jet and 
water jet impact onto a flat target produces the same 
eject a sheet angle ^o, the experimental results presented 
here reveal that this similarity obtains in the presence of 
highly dissimilar internal structures. Granular jet impact 
produces a dead zone while water jet impact creates an 
axisymmetric straining flow without a dead zone. This 
makes the similar ^o values puzzling. 

To understand this, we recall the observation from 
Cheng et al. that, if the entire jet is modeled as a sin- 
gle degree-of-freedom system with an average velocity, in 
analogy with particle impact in the dilute regime, then 
momentum and energy conservation requires cos ^o = 

I -{A- 5)(i?Tar/i^Jet)' whcu i?Tar < i^Jet HE]. The 

dimensionless constant A is the reaction force exerted 
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FIG. 3. (color online) Discrete-particle simulations using a 
frictional target reproduce the velocity field from experiments 
with roughened targets. Using a frictionless target in the sim- 
ulation eliminates the dead zone but does not significantly 
change the ejecta. (a) Good agreement between speed con- 
tours |u|/[/o from discrete particle simulation with a frictional 
target (left, solid), the continuum frictional- fluid model with 
a no-slip target (left, dashed) and experiment (right). The 
grey curves outline the dead zone, (b) Normalized scattering 
profiles from discrete-particle simulations for a frictional tar- 
get (blue points) and a frictionless target (red points) show 
that the degree of coUimation remains relatively unchanged 
regardless of the internal kinematic features. 



by the target, normalized by the incoming momentum 
flux Pin (defined as the total incoming jet momentum 
per unit time). The dimensionless constant B relates 
Pout, the momentum flux in the ejecta, via Pout/^in = 
1 — P(i?Tar/^Jet)^- Bccausc it quantifies the dissipa- 
tion rate incurred during jet impact, P is a normalized 
drag force. In this single degree-of-freedom model, the 
fact that ^0 is similar for water and granular jets means 
A — P is similar and therefore the forces experienced dur- 
ing impact are alike. This is surprising, especially in light 
of our experimental finding that the internal flow is con- 
siderably different. 

To explore the origin of the insensitivity of ^o, we first 
constructed a discrete-particle simulation to track how 
some of the many degrees of freedom present in dense 
granular impact evolve. In our scheme [23 , the particles 
are modeled as hard spheres that experience dynamic 
friction. Upon collision, the spheres lose a fraction of 
their kinetic energy. We consider two targets: (i) a fric- 
tional one where we decorate the target with a layer of 
stationary grains and (ii) a frictionless one where grains 
experience specular reflection upon collision with the tar- 




yield the following governing equations for the velocity 
field u(x, t) and the pressure field p(x, t) 



V-u = 

p(St + u- V)u = V-tT, 



(1) 

-pl + /iP7/|7| (2) 



FIG. 4. Solutions to the continuum granular-fluid model 
shows that changing the dead-zone size only affects the ejecta 
slightly, (a) As /i decreases, the normalized dead-zone height 
i^Dz/i^Tar contiuuously shrinks and vanishes at /x = 0. (b) 
The dead-zone aspect ratio Hdz/Rdz remains approximately 
constant for different dead-zone sizes, (c) The ejecta angle ^o 
changes only slightly even as the internal dead-zone shrinks 
from a signiflcant region to 0. (d) The dimensionless reaction 
force exerted by the target on the jet A barely changes over 
the entire range of /x. 



get. 

Frictional-target impact simulations reproduce the 
salient experimental results. As in the experiment, par- 
ticles leave the target in a thin sheet while an approxi- 
mately conical dead zone forms at the target. The nor- 
malized velocity contours from experiment using a rough- 
ened target and the frictional-target simulation agree 
quantitatively near the target in Fig. [3](a). 

Changing parameters, such as the coefficient of resti- 
tution or the friction between grains, produces only weak 
variations in the ejected sheet or the dead zone [24]. On 
the other hand, using frictionless targets instead of fric- 
tional ones produces a qualitative change: the dead zone 
is eliminated entirely. This dramatic change in the in- 
ternal state of the jet produces only a slight change in 
the ejecta: ^o changes from 37° for a frictional target 
to 45° for a frictionless one. The degree of collimation 
also remains similar (Fig.lsFb)). These results show that 
the granular ejecta is remarkably insensitive to whether 
the impact zone is static or freely fiowing but do indicate 
why granular jet impact behaves like water jet impact. 

To address this question, we analyze granular impact 
in the continuum limit. Since both experiment and sim- 
ulation show that the grains remain densely packed with 
a low effective temperature, we assume that the collec- 
tive motion in the jet is incompressible and isothermal. 
In addition, we assume that the deviatoric stresses obey 
generalized Coulombic friction: they lie along the shear 
direction with a magnitude equal to the pressure multi- 
plied by the dynamic- friction coefficient ja. In contrast, 
Newtonian ffuids, such as water, have deviatoric-stress 
components that are proportional to the strain rate but 
independent of pressure. 

Incompressibility and momentum conservation then 



where a is the stress tensor, I is the identity matrix, 
7 = Vu + Vu^ is the rate of strain tensor, and |7| = 
\/(7 • 7)/2. During steady impact, the a priori unknown 
jet surface satisfies the free stress condition an — 0. 
This boundary condition, together with a zero velocity 
condition at the target, completes the mathematical for- 
mulation of steady-state impact. 

Previous works [251429] modeling dense granular fiow 
as a frictional fiuid have found the best agreement with 
experiments by allowing the dynamic friction /i to de- 
pend on a ratio of two timescales: a microscopic parti- 
cle rearrange timescale determined by the local confining 
pressure, and a macroscopic timescale related to large- 
scale shear. In general, wherever the pressure and ve- 
locity fields vary, this ratio varies as well. Granular jet 
impact, however, is particularly simple because the lo- 
cal confining pressure is generated by impact alone and 
thus scales as pU^. As a result, the ratio of time-scales 
remains essentially uniform over the impact region. We 
are therefore able to reproduce the experimental mea- 
surements by choosing /i to have a single value. 

We use the open-source, time-dependent free-surface 
solver Gerris [26] [SO] [Sj to obtain steady-state solutions 
to the frictional fluid jet impact problem. The solu- 
tions reproduce the pertinent features of the experiment 
and the discrete-particle simulations. Imposing no-slip 
boundary conditions at the target, so that the tangential 
speed of the fluid is zero there, creates a conical dead- 
zone within the jet. Imposing free-slip boundary con- 
ditions, corresponding to making the target frictionless, 
eliminates the dead zone. We chose \i = 0.45 to repro- 
duce the experimental ejecta angle ^o- This produces 
velocity fields that agree exceptionally well with experi- 
ment (Fig.[3](a)). 

Fig. Ba) shows that as /i decreases, so that the 
frictional-fiuid model approaches perfect-fiuid fiow, the 
dead zone shrinks continuously and vanishes at /i = 
0, while leaving H^z/Rbz approximately constant 
(Fig. [4](b)). In Fig.[4](c), we observe that ^o only varies 
by 20° even though the internal structure has changed 
dramatically. The variation is due almost entirely to the 
reduction in the normalized drag force B. The normal- 
ized reaction force A is nearly constant in /i, regardless of 
whether impact creates a sizable dead-zone or a vanish- 
ingly small one as shown in Fig.Qd). These results sug- 
gest that the insensitivity of the ejecta to internal struc- 
ture is not due to dissipation or jet granularity; rather it 
owes its origin to the reaction force generated by perfect- 
fluid flow being insensitive to internal structure. 

To conflrm this hypothesis, we consider the impact of 




FIG. 5. (color online) Perfect-fluid calculations demonstrate 
insensitivity of impact dynamics to the target structure, (a) 
The ejecta angle ^o versus the inclination angle ^ of a coni- 
cal structure placed at the target center (Fig. Is] (b)), where 
^ < (blue) is a depression into the target and ^ > (red) is 
a conical dead-zone-like protrusion, demonstrate that even 
values considerably larger than the dead zone in our experi- 
ment show little signature in the ejecta. (b) Nearly identical 
jet surface cross-sections (dashed lines) for perfect-fluid im- 
pact against different targets (solid lines) nearly coincide. 



contrast with recent theory and simulations by Sano and 
Hayakawa [32] , our experiment clearly shows that a large 
dead-zone forms on impact. Furthermore, our simula- 
tions show that thin ejecta sheets form generically when 
the effective temperature is low and the density is high, 
regardless of whether or not an interior dead-zone is 
present. Continuum modeling shows that the dynam- 
ics are well approximated as an incompressible frictional 
fluid. Changing the dissipation produces little variation 
in the ejecta and the target reaction force, even though 
the dead zone size varies greatly. This persists to the lim- 
iting case of a perfect fluid, which is also insensitive to 
the internal structure. Thus the ejecta formed by granu- 
lar impact is robust because they are dictated by incom- 
pressibility and inertia. 
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Kadanoff-Rice fellowship, and J. E. acknowledges sup- 
port from a University of Chicago Millikan fellowship. 



a perfect fluid onto a conical structure of inclination an- 
gle 0. Positive ^ is a protrusion, mimicking the dead- 
zone geometry, and negative 6> is a depression. We match 
the cone-base radius to the experimental dead-zone ra- 
dius, 0.76i?Tar7 and vary 6 (Fig. IsFa)). The figure shows 
that ejecta sheets from three very different targets are 
nearly identical. Fig. p](b) shows that changing 6 be- 
tween —30° and +45° produces a variation in the ejecta 
angle of merely 6 percent. The experimental dead-zone 
corresponds ioO = +27°, within the range oiO presented. 

Taken together, our results demonstrate that, ab- 
sent direct measurements of internal states, a collimated 
ejecta pattern is just as plausibly produced by far-from- 
equilibrium collective motion as from a fully thermalized 
Newtonian flow. Here, collimated ejecta coexist with the 
creation of a cold dead-zone. Thus, the ejecta cannot be 
easily used to determine the internal state of the dense 
beam. This is relevant to the elliptical flow observed in 
the quark-gluon plasma at RHIC, which has been inter- 
preted as evidence for fully thermalized Newtonian flow. 

Our results are also relevant to accretion by dense gran- 
ular impact. Even without cohesive forces between the 
frictional particles, jet impact accretes an interior dead- 
zone. This is relevant to the formation of planet esimals 
which have been difficult to model [T6"-18 . Our finding 
supports the previously proposed view that the difference 
in porosity in colliding dust aggregates is more impor- 
tant than the precise strength of inter-particle cohesive 
forces [laiig. 

In conclusion, we have investigated the relationship be- 
tween the highly collimated, thin-sheet ejecta and the 
internal dynamics of a granular jet during impact. In 
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